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We theoretically investigate a scheme for atom interferometry that surpasses the standard quan- 
tum limit. A four-wave mixing scheme similar to the recent experiment performed by Pertot et al. 
[2(il ] is used to generate sub-shot noise correlations between two modes. These two modes are then 
interfered with the remaining two modes in such a way as to surpass the standard quantum limit, 
whilst utilising all of the available atoms. Our scheme can be viewed as using two correlated inter- 
ferometers. That is, the signal from each interferometer when looked at individually is classical, but 
there are correlations between the two interferometers that allow for the standard quantum limit to 
be surpassed. 



I. INTRODUCTION 

Atom interferometry [l| is a useful tool for sensitive 
measurements of gravitational fields and accelerations 
0) El > gravitational gradients [H , rotations 01 , magnetic 
fields @, the gravitational constant (G) 0, and the 
fine structure constant (a) @. Although most state 
of the art atom interferometers currently utilise laser 
cooled atoms, there may be some benefit to using Bose- 
condensed atoms. Bose-Condensed atoms provide im- 
proved visibility in configurations which require complex 
manipulation of the motional state, such as high momen- 
tum transfer beamsplitters Q . The ultimate limit of the 
sensitivity of any interferomctric device which uses un- 
corrected particles in each arm is the standard quantum 
limit A6 = -As=, where N f is the total number of de- 

tected particles [lOj. 

Recently, there has been much interest in surpassing 
this limit in atom interferometers via the use of quan- 
tum entanglement, with two recent proof of principles 
experiments demonstrating sub-shotnoise phase sensing 
[Til [l2| . These experiments exploit the s-wave scatter- 
ing to provide one-axis twisting [l3|, [l4| to generate spin 
squeezing in the two atomic modes. However, both of 
these experiments use only a small number of atoms (200- 
1200 atoms), so the absolute sensitivity of the device is 
low. It is unlikely that the number of atoms used in 
this scheme can be increased significantly, as the phase 
correlations produced become very sensitive to classical 
uncertainty in the total number of particles as the total 
number of particles is increased. 

Atomic four- wave mixing [Pol ] has long been considered 
a possible method of creating entanglement between spa- 
tially separated atomic modes (la - [23| , and the generation 
and detection of quantum correlations using this process 
has recently been demonstrated [U, [25| . Four- wave mix- 
ing has an advantage over one-axis twisting schemes as 
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the correlations produced by the interaction are number 
correlations, and are not as sensitive to the total num- 
ber of particles. However, it is so far unclear how best 
to utilise these correlations for the purpose of atom in- 
terferometry. We consider a setup similar to the recent 
experiment of Pertot et al. [2(| . Here, collisions between 
two electronic states allow for collinear four-wave mixing, 
which can be well represented by a four-mode model. A 
simple theory predicts that the two modes with low occu- 
pation number display sub-shot noise quantum correla- 
tions. However, in order to increase the sensitivity of our 
device, we devise an interferometry scheme which utilises 
all of the available atoms, such that it is possible to sur- 
pass the standard quantum limit while retaining a larger 
number of atoms. 

The structure of this paper is as follows. In section [TTl 
we will describe our proposed scheme. In section IIII1 we 
will introduce a simple four-mode model to calculate the 
level of quantum correlations generated from the four- 
wave mixing process. In section IIVI we will introduce a 
more realistic multi-mode model that allows for effects 
such as spatial inhomogeneity and multiple order scat- 
tering. In section IVl we will demonstrate how this system 
can be used to perform an interferometric phase measure- 
ment which surpasses the standard quantum limit, and 
in section I VI I we will compare our scheme to one- axis 
twisting for a large number of particles. 



II. FOUR- WAVE MIXING. 

Our scheme involves two stages: The state prepara- 
tion stage (Fig. (HJ), and the interferometry stage. 
We begin the state preparation stage with a BEC con- 
taining approximately 2 x 10 5 87 Rb atoms, tightly con- 
fined in the z and y direction, and weakly confined in 
the x direction. The x confinement is switched off be- 
fore the state preparation begins. All atoms are in the 
| a) = \F = l,mp = 1) hyperfine state. To induce the 
four-wave mixing, we couple approximately 50% of the 
atoms to the \b) = \F = 2, mj? = 0) hyperfine state via a 
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FIG. 1: (Color online) Four-wave mixing, s-wave collisions 
between |a,0) and \b, ko) scatter atoms into \a, ko) and \b, 0). 



two-photon Raman transition, which also transfers linear 
momentum hko to the atoms. In addition, we use reso- 
nant microwave coupling to transfer a small amount of 
population (approximately 1%) between \a) and |6) with 
no momentum transfer. At this stage, modes |a,0) and 
\b, ko) each contain approximately 10 5 atoms, and modes 
|6, 0) and |a, fco) each contain approximately 10 3 atoms, 
where the notation \a(b),k) = \a(b))\k), and the state 
|fc) indicates a wavepacket traveling with mean momen- 
tum hk. s-wave collisions between |a, 0) and |6, fco) atoms 
will transfer particles from |a,0) to |a, fco), and \b, fco) to 
| £>, 0) . The creation of one \a, fco) atom indicates that a 
\b, 0) atom has also been created, as well as the removal 
of one atom each from modes |a, 0) and \b, fco). This is the 
source of the quantum correlations. The 1% seed is op- 
tional, as the tight transverse confinement strongly sup- 
presses scattering into other modes. However, we found 
that the addition of a small seed improved the multimode 
dynamics of the system, as will be discussed in section lPVl 
After an amount of time tf wm , the desired level of four- 
wave mixing has been achieved, and the two momentum 
modes completely spatially separate, such that at time 
ti, \a,0) and \b, 0) are spatially overlapping, centred at 
locations xl, and |a, fco) and |6, fco) are spatially overlap- 
ping, centred at locations xr. tf wm can be adjusted by 
adjusting the longitudinal and transverse trapping fre- 
quencies, which will adjust the ratio of the characteristic 
four-wave mixing time and the time taken for the two 
momentum wave packets to spatially separate. Alterna- 
tively, a Feshbach resonance could be used to switch off 
the s-wave interactions after an amount of time tf wm . At 
this stage, the system can be thought of as two spin-i 
systems. One, formed by |a, 0) and |6, 0) centred at Xl, 
and the other formed by |a, fco) and fco) centred at xr 



III. FOUR-MODE MODEL 



In this section we will introduce a simple four-mode 
model to calculate the degree of quantum correlations 
generated from the four-wave mixing process. Beginning 
with a full multimode model describing the system, the 



U = E/ ^j(rWi(r) d 3 r 



(1) 



where j = a,b, and H — + V(y,z), and C/y 

. where is the inter- and intra- species scatter- 
ing length, and m is the mass of the atom. We can gain 
insight into the behaviour of the system by introducing 
a simple four-mode model by making the approximation 

4> a (r) « a y(x,y,z) + a ko y(x,y,z)e lkoz 

Mr) « b y(x,y,z) + b ka y(x,y,zy koz , (2) 

where ty(x,y, z) is the (normalized) ground state of the 
confining potential, before the longitudinal potential was 
switched off. Substituting this into Eq. ([!]) gives 
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where hu k — Hujo + f^", and huio is the ground state 
eigenvalue of H . The characteristic four- wave mixing 
rate is \ = Uq j \^(x, y, z)| 4 d 3 r, and we have assumed 
that U aa = Ubb = U a b, which is a good approximation for 
87 Rb. We have assumed that / \$f(x, y, z)\ 2 e lk ° x d 3 r sa 
0, which is a good approximation when the momentum 
spread of the condensate is small hko, or fco > —, where 
r x is the characteristic size of the condensate in the x 
direction. 

We can gain some insight into the behaviour of this 
system with a simple analytic model. The Heisenberg 
equations of motion for the two lowly occupied 'seed' 
modes are 



id ko = (u k[) + xN t ^ a ko + xaoh bl 



where Nt 
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t = a ao + b^ba + d\ o a ko + b\ a b kl 
number of atoms, which we will treat as a constant N t — > 
Nt, as it is a conserved quantity. We define the number 
of particles in each mode as N a L = a^ciQ, N^l = b^bo, 
NaR = al Q a koi and N bR = b\b ko . As (N aL ),{N bR ) > 
(N r), (NbL,) at t = 0, we can make the approximation 
a , b ko -> l3 ka e^ lu)k o\ where a = J {N aL ), and 
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b], b kn is the total 
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By ignoring the depletion from these 
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modes, which is acceptable for times xN t t <C 1, we find 
the solution to Eq. (g]) and Eq. flS)) is 



[28j |. The stochastic differential equations are 



a ko (t) = afe o (0) coshr - i6j(0) sinhr 



b (t) — &o(0) coshr — ia\. (0) sinhr , 



(6) 
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with r = x a oPk t and we have made the transformation 
flfco = afc e I (" fc « +JVtX )*, ~b Q = 6 oe »(-o+JV tX )t 

The quantity we are interested in is the relative number 
difference variance, which we define as 



V(N-) 



(8) 
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(Ni + Nj) 

Vi_j = 1 for two independent Glauber coherent states 
|27|, which is a good approximation of the quantum 
statistics of a coherently split BEC. Taking the initial 
quantum state of modes |a,A»o) and \b, 0) to be inde- 
pendent Glauber coherent states with mean populations 
(N aR (0)) = (N bL (0)) = N , we find that 



(N aR {t)) = (N bL (t)) = [N 



ijcosh2r-i, (9) 



and 
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Figure ^ shows (N aR (t)) and v a B„bL as a function of 
time. v a R,bL decreases from 1 exponentially with r. How- 
ever, the population grows exponentially, suggesting that 
our model is invalid for long times, and a model which 
takes into account depletion is required. 

A more complicated analysis of our system which 
includes effects such as quantum depletion and Kerr- 
dephasing due to uncertainties in the particle number can 
be achieved by using a stochastic phase-space method. 
Specifically, we use the truncated Wigner (TW) approach 
[28| . which we will now describe. The master equation 
for the system is found from Eq. ([3]) , and then converted 
into a Fokker-Plank equation (FPE) by using the Wigner 
representation. This equation can then be converted into 
a set of stochastic differential equations, which we solve 
numerically. By averaging over many trajectories with 
different noises in the initial conditions, expectation val- 
ues of the quantum field operators can be calculated. 
When converting our FPE to stochastic differential equa- 
tions, we ignore terms with third-order derivatives in the 
FPE, as these terms don't have a simple mapping to the 
stochastic differential equations, and can be assumed to 
be negligible when the field has high occupation numbers 



ia 
ia ko 

iPk 



{Lo + xNt)a + xl3* k Poa k (11) 

(u ko + X N t ) a k + xPoPk^o (12) 
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where N t = a* a + a* ko a ko + PqPo + P* ko P ko - 2 repre- 
sents the total number of particles. The subtraction of 2 
is required due to the correspondence between the com- 
plex variables in the stochastic differential equations, and 
symmetrically ordered operators. These complex vari- 
ables correspond to our original operators: 
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tial conditions for each trajectory of the evolution these 
equations was chosen such that they correspond to the 
quantum statistics of the initial state of interest. Ex- 
pectation values of symmetrically ordered combinations 
of field operators correspond to the ensemble average of 
these complex variables over a large number of trajec- 
tories. For example: |<x,| 2 = \{a}-a,j + hja}-), where the 
overline represents the mean of a large number of trajec- 
tories. 

We solve our equations of motion choosing our initial 
conditions as Glauber coherent states for each mode, with 
(N aL {0)) = (N bR (0)) = 10 5 , and (N aR (0)) = (N bL (0)) = 
10 3 . We let the system evolve for an amount of time 
tfwm- In practice, tf wm is determined by the separa- 
tion time for the two momentum components, and can 
be adjusted by manipulation of the longitudinal trapping 
frequency. Alternatively, x ca n be adjusted by manipu- 
lation of the scattering length via a Feshbach resonance, 
or adjusting the transverse trapping frequency. Figure 
([2]) shows the populations of each mode as a function 
of tf wm . The four- wave mixing creates correlations in 
particle number between modes \a, ko) and \b, 0), which 
we quantify by examining the normalized number differ- 
ence variance v aRt bL- As tf wm increases, v a Rfih initially 
decreases exponentially until the effects of quantum de- 
pletion limit the amount of squeezing. A minimum value 
of v aR tiL ~ 10~ 2 is reached at N t xtivrm ~ 10, before 
the correlations begin to decrease again. All other bi- 
nary combinations of modes display anti-correlations as 
a result of the four- wave mixing. The degree by which 
these quantum correlations can enhance an interferomet- 
ric measurement are show in Fig ([5]). 



The noise on the ini- 



IV. MULTI-MODE EFFECTS 

To investigate the validity of the approximations we 
made in the previous section, we introduce a one- 
dimensional model to investigate the effect of multi-mode 
dynamics on the four- wave mixing process. Assuming the 
dynamics in the y and z dimensions is trivial, we reduce 
Eq. (JIJ to one dimension by integrating over the y and 
z dimensions. We then derive a set of stochastic partial 
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FIG. 2: (Color online) (a) Populations of each mode as a 
function of tf wm . {No,l) was calculated using the four- mode 
TW model (blue solid line), and the ID multimode TW 
model (blue solid circles). (NbL) was calculated using a sim- 
ple model with no depletion (red "+" symbols), the four- 
mode TW model (red dashed line), and the ID multimode 
TW model (red squares). (N a R), and {Nbu) where omit- 
ted, as they are almost identical to {NbL), and (N a L) respec- 



tively, (b) Normalised number difference variance, 



for 



{**! j} ~ {aR, bL} (blue "+" symbols: Simple mode with no 
depletion; blue solid line: Four-mode TW model; blue solid 
circles: ID multimode TW model), {aL,bR}, (red dashed 
line: Four-mode TW model; red squares: ID multimode TW 
model), and {aL, bL} (green dot-dashed line: Four-mode TW 
model; green triangles: ID multimode TW model) . {aR,bR}, 
{aL, aR}, and {bL, bR} are omitted, as they are almost iden- 
tical to {aL, bL}. The multimode calculation agrees well with 
the 4-mode calculation for small values of tf wm , but doesn't 
reach the same degree of squeezing. 



differential equations for the complex functions ipi (x) cor- 
responding to the quantum operator ipi(x) via the same 
method as described in the previous section [18 1. 
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where n t {x) = (\ip a {x)\ 2 - jk^) + (\tpb{x)\ 2 - 535) is the 
total density, where Ax is the grid spacing. The subtrac- 
tion of is to compensate for the mean field of the 
vacuum, which is non-zero in the truncated Wigner ap- 
proach. U\d = -^jr is the one-dimensional interaction 

irr 

strength, and ro is the characteristic transverse width 
of the condensate. The initial condition is chosen to be 
analogous to the initial condition in the four-mode model: 

Mx,0) = VW^oix) + VW^o(x)e lkoX + 

V Ax 



Mx,o) 



'I0 3 * (a0 + V I0^ (x)e lko 



where m (x) are complex Gaussian noise functions with 
standard deviation in the real and imaginary components 
of one half and (rf n {xi)r\m{xj)) = h^i,j^n,m- I n the TW 
method, the expectation value of the density of state \j) 



At t = 0, con- 



atoms is {i)](x)$j{x)) = \ip 3 {x)\ 2 - 2E 
finement in the x-direction is switched off, whilst keep- 
ing the transverse confining potential, and the the two 
wavepackets are allowed to spatially separate. After an 
amount of time tf wm , the s-wave scattering length is 
switched to zero using a feshbach resonance, to control 
the amount of four-wave mixing that will occur. After 
70 ms, the wavepackets have completely separated, and 
we measure the number of atoms on either side of a cen- 
tral position xq. Using this measurement, we can define 
the quantities analogous to the populations of the four 
modes in the previous section as 
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Figure ^ (a) shows the populations in each of these 
modes as a function oitf wm . For early times NtXtfwm 
4 (tf wm ^ 150 /is), the multimode model mimics the sim- 
ple four-mode model. This is also the time when the vari- 
ances (Figure ^ (b)) begin to deviate from the simple 
four-mode model and v a L,bR reaches its minimum value 
v a L,bR ~ 0.1, almost a factor of ten away from the min- 
imum predicted from the simple four-mode model. The 
discrepancy is due to complications arising from mul- 
timode effects. Figure (]3]) shows the expectation value 
of the density of the atomic clouds after separation for 
tfwm = 0.1 ms (Ntxtfwm ~ 2.68) and tj wm = 0.18 
ms (Ntxtfwm ~ 4.82). In (a), the density profile of all 
wavepackets closely resembles that of the initial state, 
suggesting that Eq. ^ is a good approximation in this 
case. However, (b) shows significant deviation from these 
spatial modes. One cause of this is that the density is 
higher in the center of each cloud, so the four- wave mix- 
ing occurs at a greater rate there, increasing the deple- 
tion. It is interesting to note that discrepancies begin 
to arise in all three models at roughly the same point, 
as the quantum depletion which causes the un-depleted 
pump model to break down also causes approximation in 
Eq. © to break down, as the depletion alters the shape of 
the spatial mode. The 2nd effect is the finite momentum 
spread of the wavepacket allows for spontaneous scat- 
tering to unpopulated momentum modes, which slowly 
grow. This is the cause of the 'tails' leading out from the 
condensate in Figure (|3|) (b). Without the 1% 'seed' in 
the lowly occupied modes, we observe a greater degree of 
this spontaneous scattering. However, increasing the ini- 
tial population in the seed modes degrades the degree of 
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FIG. 3: (ipt(x)i) a (x)) (blue solid line) and (ij>l(x)i> b (x)) (red 
dotted line) after 70 ms of separation time for (a): tf wm = 
0.1ms, and (b): tf wm = 0.18ms. The expectation values 
were calculated from 1200 trajectories. The trapping frequen- 
cies of the harmonic potential at t = were {Ld x ,uj y ,uj z } — 
2-7r{5, 1000, 1000}Hz. The interaction strength was scaled by 
the cross-sectional area of the condensate: U\d = -^r, with 

ro = 0.55 fim. xq was chosen as the point x = 0.4 mm, which 
is roughly halfway between the two wavepackets. 



quantum correlations. Picking the optimum level of ini- 
tial population in these modes is a trade-off between sup- 
pressing spontaneous scattering into unpopulated modes, 
and degrading the degree of quantum correlations pro- 
duced. We did not do a comprehensive search to find the 
optimum value of this initial population, but found that 
the model behaved to our satisfaction with a 1% seed. 
This effect of spontaneous scattering is not present in a 
semiclassical model, as spontaneous scattering is forbid- 
den. 



V. INTERFEROMETRY BELOW THE 
STANDARD QUANTUM LIMIT 

The obvious method for using the state generated from 
the four-wave mixing process to perform phase measure- 
ments with sensitivity greater than the standard quan- 
tum limit would be to use |a,fco) and |6, 0) as the in- 
puts to an atom interferometer, as these are the modes 
which display two-mode squeezing. However, there are 
two drawbacks to this approach. The first is that each 
of these modes are only lowly occupied, so although they 
would display sub-shot noise sensitivity when used in the 
appropriate way for an interferometric measurement, the 
absolute sensitivity of the measurement would be low, 
as the total number of particles is small. The second 
draw-back is that these modes aren't spatially overlap- 
ping, so a complicated sequence of momentum-shaping 
pulses would be required in order to construct an atom 
interferometer out of these two wavepackets. By combin- 



ing these two lowly occupied modes with \a, 0) and \b, fco) , 
which have a much higher occupation, we may be able to 
maintain the advantage of the quantum squeezing, while 
increase the total number of particles, such that the ab- 
solute sensitivity is increased. Additionally, \a, 0) and 
1 6,0) form a spatially overlapping pair, as do |a, fco) and 
|6, fco). We can combine these four modes to make an 
interferometric measurement in the following way. The 
interferometry stage is implemented by a sequence of ^ 
microwave pulses resonant with the \a) — » \b) transition, 
which provides local coupling for the spatially overlap- 
ping wavepackets. Our scheme is described in Figure 
In total, three coupling pulses are used, with a phase shift 
on \b, 0) and |6,fco) before each coupling pulse. Phase 
shifts on the left and right wavepacket need to be con- 
trolled independently. This could be implemented by the 
AC Stark shift from independent focussed optical beams. 
The first coupling pulse is used to prepare the state such 
that it is useful for atom interferometry, with the phase 
shift before the first pulse chosen such that the popula- 
tions in \a, 0) and |6, 0), and |o, fco) and \b, fco) are approx- 
imately equal after the first beam splitter. At this stage, 
the quantity S = (N aL - N bL )- (N bR - N aR ) is squeezed. 
That is, the variance in this quantity is much less than 
for uncorrelated sources. The phase shift before the sec- 
ond coupling pulse is again chosen such that the popu- 
lations remain equal. At this point, S is anti-squeezed, 
and the relative phase will be squeezed. The physical 
process of interest causes a phase shift 4>2L between |et, 0) 
and 0), and 4>2R between |a, fco) and |6, fco) to be accu- 
mulated during the free evolution between the 2nd and 
3rd coupling pulses. After the 3rd coupling pulse, the 
the number of atoms in each mode is measured, and S is 
calculated. 

For any interferometric device, the smallest detectable 
phase shift is 



A<j>- 



d{S) 
d<j> 



(21) 



where S is the measured signal. The ultimate limit for 
sensitivity is the Heisenberg limit Acj) = jr, where N t 
is the total number of detected atoms. However, when 
using uncorrelated atoms as the input to the interferom- 
eter, the maximum sensitivity is given by the standard 
quantum limit, A<ft = ^= [101 ] . 

After the wavepackets have spatially separated, we es- 
sentially have two interferometers; one centred at x^, 
and the other centred at xr. When treating our two in- 
terferometers separately, the sensitivity of each device is 
well above the standard quantum limit. However, when 
taking our signal to be the difference in the two signals: 
S = [NaL - N b L) - (N b R - N u r), it may be possible to 
observe an enhancement in sensitivity, as we have correla- 
tions between the left an right systems. We will consider 
the case where the same applied phase 2 is applied to 
the left and right interferometers, with a constant off- 
set of 7T added to the right interferometer, ie <p2L = 4>2, 
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FIG. 4: Pulse sequence. After the two momentum modes 
have spatially separated, the interferometry stage is imple- 
mented by a sequence of three ? pulses. An appropriate phase 
shift 4>ql is applied to \b, 0} and <f)oR to |6, fco) before the 1st 
coupling pulse and <f>\L and 4>\r before the 2nd pulse to en- 
sure that the populations in each mode remain approximately 
equal after each pulse. These phase shifts may be different, 
and can be implemented using the AC Stark shift from a lo- 
calised optical beam. The phase shifts 4>2L and 4>2r between 
the 2nd and 3rd pulses is the quantity that the interferometer 
measures. The population in each mode is measured after the 
3rd beam splitter. 
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FIG. 5: (a): (N a L — Nth) (blue solid line: four-mode 
TW model; blue "+" symbols: ID multimode TW model), 
{NbR — Nan) (red dashed line: four-mode TW model; red: ID 
multimode TW model), and (S) (green dot-dashed line: four- 
mode TW model; green circles: ID multimode TW model) 
as a function of (j>2- The four- mode model has slightly more 
fringe visibility than the ID multimode model. The horizon- 
tal black dotted lines are a visual guide only, (b): ffl from 

the four-mode model (solid blue line), and the multi-mode 
model (dashed blue line) as a function of <j>2. The horizontal 
black dotted lines indicated the level expected for uncorre- 

lated atoms (log 10 ( J?f) ) ~ and squeezing below the mi- 
correlated level by a factor of 10 (log 10 ( Jn) ) = — -0- ( c ) The 
phase sensitivity Acj> from the four-mode model (blue solid 
line) and multi-mode model (blue dashed line). The mini- 
mum of Acf>yNt for the multimode model is Acj) « OAl-i/Nt, 
more than twice as sensitive than for uncorrelated atoms. 



4>2R — 7T + 4>2 ■ In this configuration, the system is set up 
to measure some quantity that effects each interferom- 
eter equally, such as a homogenious magnetic field. By 
removing the phase shift 7r from the right interferometer, 
the system is now configured to measure the difference 
in the external phase shifts, and is suited to measure 
some quantity which is different for each interferometer, 
such as detecting the difference in magnetic field at two 
points. Figure ^ shows (S), an d A0 as a function 

of <p2 after the 3rd coupling pulse for tf wm = 0.12ms. 
The minimum of A<f>y/Nt for the multimode model is 
A(j)y/Nt ~ 0.41, more than twice as sensitive than for 
uncorrelated atoms. The fringe visibility of the final in- 
terferometric measurement is ~ 0.999 for the four mode 
model, and ~ 0.935 for the multimode model. 



VI. COMPARISON WITH ONE-AXIS 
TWISTING 

An alternate scheme for performing atom interferom- 
etry with sub-shot noise sensitivity is the previously dis- 
cussed one-axis twisting scheme This scheme 
involves splitting the condensate, and allowing the non- 
linear interactions to 'shear' the relative phase of the two 
modes until the desired level of spin squeezing is reached. 
However, at this point, the squeezing is not directly us- 
able, and the quadrature squeezing needs to be rotated 
by an angle 6 such that there are reduced fluctuations 
in the relative phase. In order to achieve this, an addi- 
tional phase shift and a coupling pulse is applied, in an 
analogous method as discussed in section|V]of this article. 
Significant anti-squeezing instead of squeezing will be ob- 
served if the incorrect phase shift is applied. This phase 
shift depends critically on the total number of particles, 
as the nonlinear effect which causes the phase shearing 
also causes a relative phase shift between the two modes 
which depends on the number of particles. As the num- 
ber of particles becomes large, a small variation in the 
total number of particles is sufficient to significantly al- 
ter the relative phase of the two modes, and thus rotate 
the squeezing quadrature such that anti-squeezing is ob- 
served. 

To compare the performance of the one-axis twisting 
scheme to our four-wave mixing scheme, we analysed 
the one-axis twisting scheme via a two-mode model de- 
scribed in [l4| ■ To demonstrate sub-shotnoise sensitivity, 
one would have to find the optimum value of the phase 
shift, and then perform many shots of the experiment 
to build up quantum statistics in the signal. We set 
the total number of particles as 2 x 10 5 , and optimis- 
ing the nonlinear interaction time, coupling pulses, and 
phase shifts, such that the phase sensitivity of the device 
was A(f>y/7Jt R3 0.4. We then varied the total number 
by a small amount, whilst keeping all other parameters 
(phase shifts and coupling pulses) fixed. We found that 
by altering the total number of particles by about 1% 
was enough to decrease the sensitivity to the limit set by 



7 



shot-noise (A<j) = 77^)- It is important to note that this 
degradation of the signal is not due simply to the fluctua- 
tions in the number of atoms directly adding noise to the 
signal, as the signal (ie, difference in the number of atoms 
in each mode) is insensitive to fluctuations in the total 
number of atoms. A small shift in total number causes a 
big shift to the value of the phase shift required in order 
to rotate the squeezing to the correct quadrature. How- 
ever, when the total number of atoms was small (~ 1000), 
as in the recent experiments by [ill) E3 > a small fluctua- 
tion in total number causes a relatively small phase shift, 
and a perturbation in the number of particles of ~ 35% 
was required in order to decrease the sensitivity to the 
shot-noise limit. 

As the Hamiltonian for the fourwave mixing process 
doesn't rely on asymmetric scattering lengths between 
the two modes, the process that produces the correlation 
does not also cause a number dependent phase shift. We 
found that altering the initial number by 50% in a simple 
four-mode calculation, while keeping all other parameters 
the same as was used in Figure© (x, tf wm , <f>oR, 4>ir, 
4>2R., 4>ol, 4>il, and 4>2h), the shot-noise limit was still 
surpassed. This indicates that four- wave mixing is suf- 
ficiently robust against shot-to-shot number fluctuations 
that it can be used with a large number of particles. 



VII. CONCLUSION 

We have demonstrated that quantum correlations gen- 
erated from four-wave mixing can be utilised to perform 
atom interferometry with sensitivity greater than the 
standard quantum limit. Over longer timescales, multi- 
mode effects cause complicated evolution in the spatial 
mode of the wavepackets, which reduce the strength of 
the useful correlations and ultimately limit the amount 
of squeezing once can produce. By combining lowly occu- 
pied, yet correlated modes with the two remaining highly 
occupied modes, it is possible to surpass the shot-noise 
limit, while still using a large number of atoms, increasing 
the absolute sensitivity of the device. 
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